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ABSTRACT 

Available proper motion measurements of Milky Way (MW) satellite galaxies 
are used to calculate their orbital poles and projected uncertainties. These are 
compared to a set of recent cold dark-matter (CDM) simulations, tailored specifi- 
cally to solve the MW satellite problem. We show that the CDM satellite orbital 
poles are fully consistent with being drawn from a random distribution, while 
the MW satellite orbital poles indicate that the disc-of-satellites of the Milky 
Way is rotationally supported. Furthermore, the bootstrapping analysis of the 
spatial distribution of theoretical CDM satellites also shows that they are con- 
sistent with being randomly drawn. The theoretical CDM satellite population 
thus shows a significantly different orbital and spatial distribution than the MW 
satellites, most probably indicating that the majority of the latter are of tidal 
origin rather than being DM dominated sub-structures. A statistic is presented 
that can be used to test a possible correlation of satellite galaxy orbits with their 
spatial distribution. 



Subject headings: Galaxies: kinematics and dynamics; Galaxies: Local Group; 
Galaxies: evolution; Galaxies: haloes 
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Introduction 



The Milky Way (MW) is surrounded by about 20 known satellite galaxi es of which the 



very faint ones were found most recently by systematic search campaigns 


kVillman et al. 




2005 


Sakamoto & Haseeawa 2006: Zucker et al. 


2006a 


b; 


Belokurov et al. 2006 


. 2007; 


Walsh et al. 


2007 


) in the Sloan Digital Sky Survey (SDSS, 


York et al 


I2OOO) catalogue data. These satel- 



lite galaxies are unique laboratories to study galaxy evolution since they are close enough to 
resolve individual stars, they are massive enough to have possibly kept back gas in the past, 
fuelling star formation processes, and they are old enough to study formation processes of 
galaxies in the early Universe and their evolution over a Hubble time. Yet, the true nature 
of the dwarf-spheroidal (dSph) satellites is still unclear. 

The dSph satellite galaxies of the MW (and Andromeda) are mostly believed to be 
lumin ous, cold dark-matte r (CDM) dominated subhaloes that fell into the Milky Way's host 
halo (IWhite &: Reed Il978l ). A lot of effort has recently been made to push the limits of 
CDM simulations to make predictions of luminous CDM sub-structures within MW sized 
host-haloes. Performing hierarchical structure-formation simulations combined with algo- 
rithms describing dark-matter - luminous-matter biases, as given by semi-analytic galaxy 
formation models for example, DM sub-ha l oes are generated that a ppear to resemble dSph 
satellite galaxies (e.g. iLibeskind et al.ll2005l . 120071 : iMoore et al.ll2006l ). Given the discrepancy 
between the observed visible matter and the inferred total mass, dSph sa tellites would be 



the most dark-matter do minated objects currently known in the Universe ( iWilkinson et al. 



2OO2I : iKleyna et al.l 20031). Their inferred dark matter profiles, however, are incompatible 
with CDM theory (IWilkinson et al.ll2006l ). Recent observational constraints question the 



validity of the CDM model in these systems, and imply that perhaps t he warm dark-matte r 
(WDM) model may need to be considered as a viable alternative ( IGilmore et al.l 120071 ). 
This interpretation relies on solving the Jeans equation which implies that the dSph satel- 
lites have cores extending to about 100 pc, wher eas CDM haloes are known to be cuspy 
with a p oc r~^, p ^ 1 - 1.5, density profile (IDiemand et al.l l2005l ). Independent evi- 
dence for the existence of cores come from observations of sub-structure within Ursa Minor 



(iKlevna et al.l 120031) and th e placement of globular clusters in Fornax fiGoerdt et al.l 12006 
Sanchez-Salcedo et al.ll2006l ). 



An alternative theory for the origin of dSph satellite galaxies of the M W is that they are 
of tidal origin: given the s atellites fall on one o r two great circles on the sky (iKunkel fc Demers 
19761 : lLynden-Belllll976l ). iLynden-Belll (119831 ) first suggested that the satellite galaxies of the 
MW might be torn off a common progenitor galaxy. With the advent of large N-body CDM 
simulations this idea was often overlooked, since in CDM simulations dSph galaxies are 
naturally identified with accreted cosmological sub-structures, even though not all of their 
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properties could be fully understood (see above). iKroupa et al.l (120051 ) highlighted this issue 
again by examining the anisotropic spatial distribution of the MW dSphs in the context of 
CDM, pointing out tha t their distribu tion can be best represented by a plane highly inclined 
to the MW stellar disc. iKroupa et al.l deemed this to be incompatible with dark-matter dom- 
inated satellites, pointing out instead that the plane suggests a causal connection between 
the satellites. They emphasised that a correlation between the satellites would naturally 
arise if they were of tidal origin, i.e. tidal-dwarf galaxies, originating from an interaction 
with the young MW. 



Tidal-dwarf galaxies (TDG. iMirabel et al.lll992l ) should be found frequently: on the or- 
der of ten or more TDGs have been observed to f orm in single galax y interaction event 



in the local Universe (e.g. IWeilbacher et al.l l2003l : IWalter et al.l |2006| ). The production 



rate in t he early Universe is expected to be even higher given the larger gas fraction in 
galaxies. lOkazaki fc Taniguchil (120001 ) argued that if only 1-2 long-lived TDGs form per 
encounter, the number of dwarf ellipticals (dEs) in very different environments can be ex- 
pl ained as being tidal-dwarfs. T his number of TDGs per encounter is deduced independently 
by iDelgado-Donate et al.l (l2003l ) from direct observations strengthening the case for TDGs 
being a major contributor to the dwarf galaxy population. The formation of TDGs has also 
been seen in high resolution simulations of gas-rich encounters of galaxies (iBournaud &: Due 
20061 : IWetzstein et al.l 120071 ). although very little is known about their fate. It has, however, 
been sho wn that dark-m atter free dwarf galaxies can survive for several Gyr in a MW like 
host halo (lKroupalll997l ): these computations have shown the existence of long-lived solutions 
in phase-space that appe ar dark- matter dominated but are instead d erived from DM free 
satellite galaxies (see also iFleck fc KuhnI l2003l : iMetz fc Kroupal 120071 ) . The essential point 
to note is that TDG-formation is an inherent part of any cosmological structure formation 
theory, and sufficient stellar-dynamical and observational work already exists proving their 
long-term survival. So the natural question to ponder is: where are they in the Local Group? 

To constrain the possible origin of dSph satellite galaxies one can study internal processes 
like radial density- and velocity-dispersion profiles, star-formation histories, or the chemical 
enrichment of stars. Supplementary to this, their spatial distribution and orbital properties 
can be investigated. The spatial distribution of the MW satellites has been shown to be 
highly anisotropic: the satellite galaxies are arr anged in a disc-like structure, lying almost 



perpendicular to the Milky Way's galacti c disc (iKunkel fc Demerd Il976l : iLynden-Belll Il976 



19821 : iHartwickl I2OOOI : IKroupa et al. 



2005: Metz et al 



200 



3)^ This disc-like alignment, the 



disc-of-satellites (DoS), is statistically very sign ificant and incompatible with being ra ndomly 
drawn from a distribution of CDM sub-haloes ( Kroupa et al. 20051: Metz et al. 2007). Based 
on th e spatial alignment, iKunkel Sz Demers Jl976h . iLynden-Belll Jl976h and iLynden-Belll 
( 119821 ) postulated that some of the dSph galaxies, along with several globular clusters (GC), 
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belong to one or two streams of co-moving objects. iLynden-Belll (119831 ) first suggested that 
these co-moving dSph galax ies might be tidally torn off a com mon progenitor galaxy. In- 
corporating radial velocities, iLynden-Bell &: Lynden-Belll (119951 ) further constrained possible 
streams of dSphs and GCs and used t he great circ l es of possible locations of angular mo- 
menta to identify possible correlations. iPalma et al.l (120021 ) extended the analysis of these so 
called "polar-paths" by using available proper-motion measurements to confine the location 
of directions of the angular momenta of satellite galaxies to arc-segments of the polar-paths. 

Meanwhile further proper motion measurem ents, especial l y for more d i stant dSph galax- 
ies, have become available. In a series of papers iPiatek et al.l (l2003l . |2005| . |2006| . 120071 ) mea- 
sured the proper motions of the dSph galaxies Carina, Ursa Minor , Sculptor an d Fornax, 



respectivel y, with the Hu bble Space Telescope (HST). Here we follow iPalma et al.l : we mcor- 
porate the iPiatek et al.l measurements to construct new arc-segments of polar-paths for the 
MW satellite g alaxies in ^2.11 Th e analysis is concentrated on the 11 "classical" (brightest) 
MW satellites (IMetz et al.l 120071 ). because theoretical (dark matter) modelling of satellite 
populations has been focused on this sample only. It is also for this sample t hat velocity 
information are available. In §2.21 we analyse, using the techniques developed in lMetz et al. 



(120071 ). the spatial distribution of a set of simulated luminous CDM haloes published with 
the claim to be a solution to the disc-of-satellite problem. The distribution of the angular 
momenta of the CDM satellite galaxies are calculated. In §2.3l the distribution of the orbital 
poles for the classical MW and the CDM simulations are statistically analysed. The results 
are discussed in O and we summarise our conclusions in ® 



2. Analysis of orbital poles 
2.1. The Milky Way satellite galaxies 



Available proper motion measurements for the satellite galaxies of the MW were col- 
lected from the literature, and are given in Table [H We do not give a complete list of all 
available measurements for the LMC and SMC here, but only the most recent ones. The 
directions of the angular momenta are not affected by this choice. Multiple measurements of 
the dSph satellite galaxies are presented where available, however each individual measure- 
ment is less precise than the most recent LMC/SMC values. We also indicate whether proper 
motions were determined using ground based techniques (typically based on old photographic 
plates with long basehnes ) or using the Hub ble Space Telescope (HST). The published proper 
motion value for Carina (IPiatek et al.ll2003l ) was derived without the advance d charge trans- 
fer inefficiency (CTI) correction for the STIS camera (IBristow et al.ll2005l ). A re-derived 
value is given in the table. For the satellites Ursa Minor, Sculptor, and Fornax a weighted 
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mean proper motion using all available measurements was calculated, taking into account 
all independent single measurements with HST. These mean proper motion values were used 
in the following analysis. We note that all measurements for a given satellite are consistent 
with each other within 3 a. 

Heliocentric space velocities were calculated from the proper motions and radial veloc- 
ities from the literature, and corrected for the circular motion of the local standard of rest 
(LSR) about t he Galactic Centre using t he lAU recommended value of Vq = 220 kms~^ cir- 
cular velocity (IKerr fc Lynden-Belllll986l ). as well as for the pecul iar motion of the S u n wit h 
respect to the LSR using the velocity components as given by iDehnen &: Binneyl ( 119981 ) . 
Positions of the satellite ga laxies were transformed to a galactocentric coordinate system 
(table 1 in lMetz et al.ll2007l ). Having galactocentric positions and velocities the directions of 
the angular momenta or "orbital poles" of the satellite galaxies of the MW were calculated. 

The given uncertainties of the proper motion measurements were used to calculate un- 
certainties of the derived angular momenta. We did not, however, incorporate uncertainties 
of any other measured quantity: distance and position of satellite galaxies, their radial ve- 
locities, the distance of the Sun from the Galactic centre, the circular velocity of the LSR, 
and the peculiar motion with respect to the LSR. All these uncertainties are typically neg- 
ligible since they do not contribute much to the total uncertainty of the directions of the 
angular momenta (see Appendix [X]). Finally, the uncertainties in the direction of the angu- 
lar momenta were projected onto a unit sphere to derive the arc-segments rep resenting the 



proje cted la uncertainties of the directions of the angular momenta (see also iPalma et al. 
2OO2I . and figure 3 therein). 



The orbital poles of the MW's satellites and their corresponding arc-segments are shown 
in Fig. [T] by the light filled symbols and solid lines, respectively. The dashed loops show re- 
gions with 15° and 30° distance from the direction of the fitted pole of the DoS f or the eleven 
classical MW satellite galaxies (/mw = 

157?3, 6mw = -12?7: iMetz et al.ll2007f ). From this 
figure one can directly identify two satellite galaxies which can not belong to a possible com- 
m on stream of s atellit e galaxies: Sagittarius and Sculptor, both having also been excluded 
by lPalma et al. J2002h as possible stream members. Sagittarius is very close to the Galactic 
disc and may have precessed sufficiently in the non-spherical potential induced by the Galac- 
tic disc, thus being far away from the direction of its initial orbit. Sculptor, however, is more 
distant and unlikely to be affected much by precession. Compared to the other satellites 
Sculptor is on a counter-rotating but nearly co-planar orbit about the MW. The star symbol 
in Fig. [T] marks the mean spherical direction, I, of the orbital poles of the possible stream 
members: LMC, SMC, UMi, Dra, Car, and For, (Zmw = 177°0, 6mw = -9.°4). The small star 
marks the position of the mean spherical direction if the LMC/SMC centre-of-mass is used 
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instead of both galaxies individually (but see iKallivayalil et al.ll2006bl ). The soli d loop gives 



the r egion with a distance corresponding to the spherical standard distance (see iMetz et al. 



20071 for its definition), Agph = 35.°4, from the mean spherical direction of the sample. The 



coincidence of I and the normal of the DoS provides strong evidence for rotational support. 

The different proper motion measurements from the literature given in Table [1] some- 
times differ by more than la. This is an inevitable weakness of currently available data. 
Most orbital poles are however not overly affected. Sagittarius is on a perpendicular orbit 
to the DoS for both /i^ values given in the Table. Likewise, Sculptor is on a counter-rotating 
orbit for both measurements, ground based and HST, respectively. Only UMi's orbital pole 
changes when incorporating the independent data, albeit then with very large projected 
uncertainties. If we take the weighted values for the two HST fields only, its pole-distance 
is 87.3° ^^^'^o, compared to 37.1° ^'^'^o for the mean weighted proper motion of all mea- 
surements. Note that the arc-segments are projected uncertainties. While the given proper 
motion uncertainty for Car is about three times smaller compared to Dra, the projected 
uncertainty is actually two times larger. Therefore, a large proper motion uncertainty does 
not necessarily result in a large uncertainty of the location of the orbital pole. 



2.2. CDM simulations 



In the recent search for an explanation of the DoS in terms of CDM substructures 
many teams have proposed solutions with varying degrees of means. One of the most recent 



works claims to be able to arrive at the DoS for all MW-like host haloes ( Libeskind et al. 



20051 ). For comparison to the Milky Way satellites ( §2.ip this set of simulated systems of 
luminous satellite haloes in a MW-type dark-matter host halo is analysed. The simulations 
and galaxy formation models are described in that work and references therein. In short, 
selected host haloes in a large cosmological dark-matter simulation were re-computed with 
high-re solution, and ga laxy formation was treated by a se mi-analytical galaxy formation 



model JCole et alijgOOol. as exten ded by lBenson et al.ll2002l ). Our data differs slightly from 



the data used in lLibeskind et al.l due to an improved semi-analytic galaxy formation model. 
We analyse a set of the eleven DM sub-haloes containing the most luminous galaxies in each 
of the six simulated MW- type host haloes. This is a somewhat artificial cut and was chosen 
by iLibeskind et al.l (120051 ) to match the number of the "classical satellite galaxies" of the 
MW. The only information used here for the analysis are the final positions and motions 
of the sub-haloes within the host halo. It is, however, neither sure that the host galaxy is 
a disc galaxy, nor how its disc would be oriented with respect to the host halo. To be a 
system comparable to the Milky Way (or Andromeda) both attributes have to be postulated 
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to match with those of the MW. 



In a first step the bootstrapping analysis as described in iMetz et aLl (120071 ) was apphed 
to the data, to compare the properties of the spatial distributions. The resulting shape- and 
strength-parameter are listed in Table [2l None of the simulated sets has as h igh values as the 



Milky Way data. Using the same nomenclature as in lLibeskind et al.l (120051 ). the simulation 
gh7 has a slightly higher shape parameter than the MW data, but the strength parameter 
is significantly lower, which means that the distribution of directions of normals fitted to 
the bootstrapped data is significantly more widespread than found for the MW. gh2 has a 
similarly strength parameter but a significantly lower shape parameter; the bootstrapped 
normals are less circular or show a multi-modal distribution. The other two data-sets with 
clustered distributions, gh6 and ghlO, are lower in both parameters, while ghl and gh3 have 
significantly lower shape parameters. Thus none of the simulated MW galaxy haloes appear 
to reproduce both the shape and strength parameters of the MW's DoS. At a significance 
level of 1 — a = 99%, all are compatible with the hypothesis to be drawn randomly from a 
spherical symmetric parent distribution of dark matter sub-haloes. 

To add further insights we compare the distributions of the orbital poles of the satellites. 
Figure [2] shows the directions of the angular momenta of all six simulations. In contrast to 
Fig. [1] for the MW data, the coordinate system can be arbitrarily chosen, since no galactic 
disc is comprised in the simulated data which would naturally define the coordinate system. 
Therefore a coordinate system has been chosen such that the axes of the coordinate system 
are ahgned with the formally fitted plane of the spatial distribution of the simulated satellite 
galaxies. The projection is chosen such that (i) it is centred on the normal, n, of the fitted 
disc - which would be the mean rotational axis if the fitted plane were rotationally supported 
- and (ii) it was chosen such that the plot is centred on the hemisphere that contains the 
mean spherical direction of the orbital poles Z, which is marked by a star symbol. So the 
plots are centred on the direction 

rii = sgn(n ■ I) n (1) 

where sgn(x) = x {xl^-^. 

The choice of the projection for Fig.[2]makes the plots also easily comparable with Fig.[T) 
If the fitted plane is rotationally supported, all or at least most angular momenta should 
lie close to the pole of the fitted plane. For comparison the same regions as in Fig. [1] are 
marked by the dashed loops. A strong clustering in or close to this region is not present in 
the simulations, and the spherical mean direction, even for selected sub-samples, is nowhere 
close to the plane's pole. The agglomeration of five orbital poles seen for the simulation ghl 
turned out to belong to a spatially distinct sub-group of sub-haloes at about 200 kpc from 
the hosts centre that have an approximately common stream motion. This motion may be 



- 8 - 



inte rpreted as an infa lling; or bypassing group of sub-haloes that has not yet disbanded (see 
also lLi fc Helmill2007l ). The remaining seven sub-haloes are all found very close to the centre 
of the host halo in a more or less spherical distribution and the directions of their angular 
momenta are widely spread. 



2.3. Pole-distance distribution 

For the MW, as well as for the simulated data, a reference direction n; (Eq. [1]) is given. 
This would be the approximate orbital pole under the hypothesis that the disc-like structure 
is rotationally supported, a hypothesis borne out for the MW as we have seen. We construct 
a pole-distance distribution (FDD) by measuring the angular distances of the orbital poles 
from this reference direction. // the disc is rotationally supported a significant excess of 
small pole- distances over a random distribution is expected. 

In Fig. [3] the cumulative pole-distance distribution for the measured angular momenta 
of the MW is plotted. For the simulations we build the sum of all individual pole-distance 
distributions. The cumulative pole-distance function expected for a random distribution, 
-Dpdd(0) = |(1 — COS0), < < TT, is shown by the dotted line. This is, however, not the 
true parent random distribution for a measured FDD, since the measured FDD is biased: 
we have chosen that direction of the normal of the fitted plane as the reference pole that 
is the presumable rotation-pole ( §2.H §2.21) . Therefore a biased FDD was constructed: 
random directions are created and the mean spherical direction r of these is determined. As 
the reference directions the x, and z axes are chosen. The pole-distance is measured from 
+x if the x-component of r is positive, from —x otherwise. The same is done for the y and 
z components. The resulting cumulative biased random FDD is shown in Fig. [3] by solid 
curves for N = 11 (CDM) and N = 8 (MW), respectively. The cumulative biased random 
pole-distance distribution, created by the Monte Carlo method, can be fitted by 

^pdd(0, AT) = 1 (1 - cos 0) + ^ sin^ (2) 
^ V A* 



for A^ > 5, < < TT. Only in the limiting case limTv^oo Eq- (ISD goes over into -Dpdd(0) = 
|(1 — COS0), which is the purely random, unbiased distribution. 

The FDDs were tested against the bia sed random po le-distance distribution using a 
Kuiper statistic V and a Watson statistic U'^ ( lStephenslll974j ) . Both statistics were originally 



used for distributions on a circle and are thus useable for the FDD. For both statistics and 
at a significance level 1 — a = 99% we can not exclude that the FDD for each individual 
CDM simulated sub-halo system is drawn from a random biased FDD (A^ = 11), nor can 
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we exclude that the co-added PDD for all simulated systems is drawn from a random biased 
PDD. The null-hypothesis can also not be excluded for the Milky Way (A^ = 8), but given 
the very large uncertainties this is not a very strong restriction. Within the uncertainties of 
the measurements we find the directions of the angular momenta of Dra and Car close to the 
spherical mean direction (< Agph). In that case we can exclude at a very high significance 
level 1 — a = 99% that the MW PDD is drawn from a random biased PDD. 

Interestingly, the CDM simulated systems show an excess of pole-distances at ~ 90°, 
which implies that the fitted disc-like structure tends to disperse. In contrast, we do find 
an over-density of small pole distances for the MW, suggesting that its DoS is rotationally 
supported. Note also that we did not choose a reference direction that minimises the pole 
distances for the Milky Way. If we had chosen the spherical mean direction, the over-density 
of small pole distances would have been more extreme. 



3. Discussion 

The bootstrapping analysis of the simulated sub-systems shows them to be compatible 
with being drawn randomly fr om a spherical iso tropic distribution, whereas this is strongly 



excluded for the MW satellites flMetz et al.ll2007l ). Moreover, the analysis of the pole-distance 



distributions (previous Section) indicates that none of the theoretical CDM sub-structure 
systems can be rotationally supported, but the orbital poles of all simulated systems are 
consistent with being drawn from a random distribution. Instead, the claimed "disc" of 
simulated sub-haloes even tends to disperse rather than being supported by rotation. The 
CDM subhalo "disc" is therefore a pressure supported flatte ned tri-axial sph e roid rather 



than a rotationally supported disc. This finding is supported by lLibeskind et al.l (120071 ) . who 
showed that the mean angular momentum of a different set of simulated satellites, using a 
smooth particle hydrodynamic code, is not co-aligned with the normal of a fitted plane to 
these sateUites. Instead it tends to be perpendicular to the long axis of t he fitted triaxial 



ellipso id, implying that the discs are tumbling around. Note, however, that iLibeskind et al. 



(120071 ) were only able to resolve galaxies brighter than My ~ — 12 which complicates direct 
comparisons with the MW satellite populations, since most MW dSph are much fainter than 
that limit. 

In contrast, the Milky Way pole-distance distribution, which is still affected by large 
measurement uncertainties, shows an excess of small pole-distances implying the DoS to be 
rotationally supported. Within the uncertainties, the two possible DoS satellites with the 
most deviating directions of the angular momenta, Draco and Carina, are consistent with 
being found close to the suggested orbital pole of the disc-of-satellites (< Agph, Fig. [T]). That 
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case would be highly inconsistent with a random sample (> 99%). In order to be able to 
draw strong conclusions based on the PDD alone we need to await more, and more-accurate 
proper motion measurements. For two satellites, Dra and Car, only a single measurement 
is available and their projected uncertainties are very large, but measurement campaigns 
with HST are in progress (Pryor, priv. comm.). With upcoming new data, also from the 
satellite missions Gaia and SIM, the statistics presented here can be used to further test the 
distribution of orbital poles. 

Another fact subs tantiates the claim of a rotationally supported DoS: the unweighted 
disc-fitting algorithm (IMetz et al.l 120071 ) is most strongly affected by the spatial location 
of the outermost satellite galaxies. We do find, however, a clear alignment of the angular 
momenta of the innermost satellite galaxies with this fitted disc, most strongly for the LMC 
and SMC. The only exception is Sagittarius whose orbit might already significantly deviate 
from its initial orbit. It would be rather curious to find a disc-like structure where only the 
inner satellites orbit within the disc, while the outer ones which mostly define the disc don't. 
Interestingly, most of the newly discovered dSph also lie close to the DoS (Metz et al., in 
preparation). It is the combination of these findings, t ogether with t he results showing tha t 
dSph can be resembled by evolved DM free galaxies (IKroupal 119971 : iMetz fc Kroupal 120071 ). 
that strengthen our notion that the dSph satellite galaxies may not originate from CDM 
sub-structures. 

Our finding for the pole-distance distribution of the Milky Way satellite galaxies poses 
another challenge to explain them as primordial cosmological substructures, adding to the 
bunch of already existing issues on galaxy-sized scales. The cusped profiles found in simu- 
lations for CDM dor ninated haloes are in disagreement with the density profiles inferred for 
MW dSph satellites (iGilmore et al.ll2007l'). but these are belie v ed to be the most d ark-matter 
dominated objects currently known (jWilkinson et al.l 120021 : iKleyna et al.l l2003l : Dabring- 
hausen, Hilker, & Kroupa, 2 008, submitted). No p hysical mechanism has been found to be 
able to evolve cusps to cores (iGnedin fc Zhadl2002l ). The missing satellite problem has been 
addressed by many groups typically explaining the difference in the number of dark-to-light 
sub-haloes by invoking a suppression of star-formation in low mass sub-haloes. Within the 
CDM paradigm, the suppression of star form ation is usually erapirica lly explained either by 



feedback processes that heat the halo gas (e.g. iKauffmann et al.lll993l ). or the e poch of reion- 



izatio n which inhibits the formation of small satellites not already seeded (e.g. iBenson et al. 
2OO2I ) . Some recent simulations were able t o reasonably reproduce the numbers and luminosi- 
ties of the bright end of the dwarf galaxies (ILibeskind et al.ll2007l ). But a complicating factor 
is that in the literature, depending on which parameters are chosen to match with, different 
studies use different criteria i n the pairing of su b-haloes with satellite galaxies. These are: 



the most massive sub-haloes (IStoehr et al.ll2002l ). the most massive progenitors of luminous 
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sub-h aloes (iLibeskind et al.ll2005l ). alternatively the earli est forming haloes (IStrigari et al. 
20071 ) , low- mass systems that were a ccreted very recently (ISales et al.ll2007b[). or sub-haloes 
that were grouped before accretion (ILi &: Helmil 120071 : iLake fc D'Onghiall2008l ). 



models by 


Besla e 


; al. 


(2007) 


( Kallivavalil et al. 


2006a 


y)t 



3) they concluded that if a ACDM-motivated model for the Milky 



Way halo is assumed, then the LMC/SMC are likely on their first passage about the MW. 
However, our results show that the orbital poles of both are correlated with the spatial 



Lvnden-Bell 


1982: 


Lvnden-Bell & Lvnden-Bell 


1995: 


Palma et al. 



2OO2I ) . Indeed, the likelihood that an uncorrelated first-time incoming satellite has an orbital 
angular momentum pointing 9° away or closer to the mean direction of the orbital poles of 
the ancient satellites (UMi, For, Dra, and Car, Imw = 179.°0, 6mw = — 4.°6) as the LMC 
does, amounts to 0.6%. It is therefore very unlikely to find such a strong correlation for 
completely unrelated satellite galaxies which they would be, if LMC/SMC are on their first 
fiy-by. Either the LMC/SMC proper-motion values are too high, the halo model is incorrect, 
or an alter native gravit ational theory has to be invoked. Among contesters for such theorie s 
are MOG jMoffatlboOsl ) and MOND/TeVeS jBekenstein fc Milgrom|[l984l : lBekensteinll2004h . 
In modified Newtonian dynamics (MOND) no object can leave the gravitational potential 
unless it is attracted by a second mass. In this case the LMC/SMC system might be bound 
to the MW, even though they seem to have a too high velocity. Interestingly, this theory 
also accounts for the extra dark-matter needed to explain observed rotation curves of TDGs, 



wher e as in Newtonian dy namics an additional mass component is needed (IBournaud et al. 



20071 : iGentile et al.l 120071 ). And, this theory would also accoun t for the ob served constant 
velocity dispersion seen among dSph satellites, a ~ 10 kms~^ ( lLokasll200ll ). whereby tidal 
and projection effects would play an important role just as in Newtonian dynamics. 

An altern ative model to expla i n the small scale phen omena mentioned above is t he tidal- 
dwar f theory (IMirabel et al.lll992l : iDuc fc Mirabellll994j ). The formation of TDGs (jZwicky 
19561 ) is observed to comraonly occur in gas-rich tidal tails thrown out of in teracting galax- 



Hunsberger et al 



1996 



Due fc Mirabel 



les (e.g. 

in simulations a s well (IBournaud fc Due 



2006 



19981: Istraughn et al. 



20061 ). and is evident 



Wetzstein et al.l 120071 ) . In contrast to the 



DoS of the MW, iBournaud fc Dud J2006h found long-lived TDGs preferentially on prograde, 
co-planar orbits, but they only considered TDGs more massive than 10^ M©, compared to 
some 10^ - 10^ Mq for the least massive dSph satellites of the MW. These models concen- 
trated on low redshift spiral progenitors. Theoretical work has also shown that dark-matter 
free galaxies can survive for a long time, and can account for most of the observed internal 
properties (luminosities, morphology, mass-to-light ratios) of dSph satellite galaxies, even 
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in Newtonian dynamics (iKroupal 119971: iMetz fc Kroupa 120071 ). Furthermore, the chemical 
evolution of TDGs has been studied (IRecchi et al.l 120071 ). suggesting that long lasting star- 
formation episodes can be understood in dark-matter free galaxies. Once TDGs have formed, 
their chemical evolution is decoupl ed from each other, which can lead to very different inter- 
nal evolutions (jCalura et al.ll2008l ). TDGs formed in one encounter, i.e. causally-connected 
TDGs, naturally align within the orbital plane of the interaction, building a rotationally 
supported disc-like structure. Their distribution may broaden due to precession of orbits 
and interactions, which a lso rnake k i ck-out scenarios of satellite galaxies due to three-body 
interaction events likely (jzhao 1998 : Sales et al. 2007a). In this theory, the only exception 
for the Milky Way satellite galaxies is Sculptor on a counter-rotating orbit, which might 
result if the tidal tails return to the host phase-shifted. More work is certainly needed to 
learn about the evolution of (old) TDGs in total: from their birth in a gas-rich early Universe 
to the possible remnants we see today. The production of TDGs must be an essential part 
of any cosmological theory - independent of how old TDGs end up: completely destroyed 
or as observable galaxies. 



4. Conclusions 

In this work we have highlighted subtle but convincing arguments that the sub-haloes 
produced in CDM simulations tailored to address the missing satellite problem can not fully 
describe the entire population of dSph orbiting the MW. Although our sample is limited 
to just six simulated galaxy haloes, we nevertheless have demonstrated an incompatibility 
between the spatial distribution of the MW satellites and the most recent CDM simulations 
that claimed to resemble the MW satellite galaxies. Additionally, the distribution of the 
orbital poles of the MW satellite galaxies also shows evidence that they differ from the CDM 
simulated systems. For the simulated haloes the orbital poles of the CDM simulations are 
consistent with being drawn from a random, isotropic distribution. In contrast, the orbital 
poles of the Milky Way show an excess of small pole-distances, but are affected by large 
measurement uncertainties. In fact, the total MW satellite angular momentum is co-aligned 
with the normal to the DoS-plane implying the DoS to be a rotationally-supported structure 
- a feature notably absent in all six simulated haloes. Such a rotationally supported disc is 
a natural consequence of their formation if the MW dSphs are of tidal origin. Indeed, the 
production of TDGs is a fundamental process in any hierarchical structure formation theory 
such that members of this class of objects must be present in the local volume. 

Clearly, the problem of the true origin of the MW satellites is of utmost importance 
for cosmological theory and further tests of models - CDM and tidal-dwarf - are required. 



- 13 - 



Here we used the currently available proper motion data for the Milky Way satellites, and 
a limited number of CDM simulated haloes. For further tests of models, and also when 
new proper motion measurements become available, a fitting function for the biased random 
cumulative pole-distance distribution is given (Eq. [2]). 

We thank an anonymous referee for her/his useful comments that helped to improve the 
paper. We are grateful to Slawomir Piatek for providing us with the re-derived proper-motion 
values of Carina. 



A. On the angular momenta uncertainties 

To calculate the projected arc-segments of the uncertainties of the direction of the 
(specific) angular momenta, I = r x v, in Section [2?T| we considered the proper motion un- 
certainties only, as these contribute to most to the uncertainties. Here we provide the estima- 



tion w hy other uncertainties are negligible. Taking equation (2) from iJohnson fc Soderblom 



( ll987n . and substituting parallaxes with distances in parsec, we get for the uncertainties of 



the velocity components in kms ^ in Galactic coordinates: 

< \ / bi2 bi3 

[{r a^^,y + {fia*(rry] \ + 2fi^,fxsk^a^ 622^23 | • (Al) 

[{r (Tf^.y + if^SCTr)^] I \&32&33 

/io*, /i5 are the proper motion components in right ascension and declination, respectively, 
given in arcsec yr~^, r is the distance in pc, and Vr is the radial velocity in kms~^. The 
constant k = 4.74057, and th e elements bij are the components of the matrix B as given by 



Johnson fc Soderblom! (119871 ). 5 is a rotation matrix, performing the transformation from 



radial velocities and proper motions to velocity components in Galactic coordinates. C is a 
matrix where each element is the square of the individual elements of the matrix B: Cij = bfj. 

The second and third component in the first term of the velocity uncertainty vector 
in Eqn. flAip are given by: (/c^ [{ra^Y + {fiarY]). With ~ O.l/i and r ^ ar, typical 
values are ar ~ O.Olr, we get A;^/i^(0.01r^ + O.OOOlr^) ^ k^af^r"^, which results in a velocity 
uncertainty of the order of ~ 25 kms^^ for each component at a distance of 50kpc. For 
the second term, with fi^* ~ fJ'S, it follows 2A;^yU^(T^ ^ k'^a'^r'^. B is a rotation matrix, with 
bij bik < I, so it is reasonable to neglect the distance uncertainties. 

Radial velocity uncertainties are typically of the order of 5 kms~^, compared to ~ 
25 kms~^ for the two transverse velocity components as estimated above. The absolute value 
of the velocity uncertainty is then estimated to be of the order 2- (25 kms~^)^ + (5 kms~^)^ ~ 
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2 ■ (25 kms~^)^. Moreover, as we calculate the specific angular momentum vector given by 
I = r X V, the contribution to the vector v from the radial velocity is almost parallel to the 
vector r. The angle is at maximum of the order arctan(8.5 kpc/50 kpc) ~ 10°, so the radial 
velocity uncertainty contributes only little to the total velocity uncertainty. 

For the distance to the Galactic Centre, we use the lAU recommended value 8.5 kpc, 
while some authors prefer a value of 8 kpc, so we assume the uncertainty to be 0.5 kpc. This 
uncertainty influences only the first component of the position vector r in Galacocentric 
coordinates. The variance of the second component of the specific angular momentum vector 
is given by af^ = (rfcx^^ + vfa'^^) + {rlcr'^^ + The contribution due to the velocity 

uncertainties is of the order 2 ■ (30kpc)^ (25 kms~^)^ = 1 125 000 (kpc kms~^)^, while the 
distance uncertainty to the Galactic Centre contributes with only about one thousandth of 
this: (100 kms^^)^ (0.5 kpc)^ = 2 500 (kpc km s~^)^. Furthermore, the distance uncertainties 
to a satellite at 50 kpc is of the order 2 kpc, so for each component of the order of 1^2^/3 kpc ^ 
1.15 kpc. These do also not contribute significantly to the total uncertainty. Finally, the 
transverse position uncertainties in a and S are much smaller than the distance uncertainties. 
At a typical distance of a satellite galaxy of the order of 50 kpc, a position uncertainty of 1' 
transforms into a spatial uncertainty of ~ 15 pc, which is two orders of magnitudes smaller 
than a typical distance uncertainty of 2 kpc. 

For the circular motion of the LSR about the Galactic Centre we use the lAU recom- 
mended value of 220 kms~^. The uncertainty of this value is of the order 20 kms~^, so this 
might have a significant influence on the location of pole of the angular momentum vector. 
We checked that the results do not change significantly by once using 200 kms~^ and once 
240 kms~^ for the circular mo tion. The uncertainties of the peculiar motion of the Sun with 



respect to the LSR as given by lDehnen fc Binneyl (119981 ) are of the order of 0.5 km s ^, much 



smaller than the other contributions and can be neglected. 

Sagittarius is the closest satellite considered here. With a distance of 24 kpc from the 
Sun, the estimations given above might not be valid anymore. But its geometric location is 
very special, as it is located almost directly behind the Galactic Centre, and r\\Vr. Moreover, 
Sgr is anyway not co-rotating with the DoS because of its spatial location with respect to 
the DoS. 
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Table 1: Measured proper motions of the satellite galaxies of the Milky Way used for the 
analysis. In the second column the galacto-centric distance of the satellite galaxies to the 
MW is given. The fifth column indicates whether values were measured using ground based 
telescopes (GND) or the HST. Weighted mean values (MEAN) are indicated as well. 





fkpc) 


( mas yr^^) 


( mas yr~^) 






Sgr 


16 


-2.65 ±0.08 


-0.88 ±0.08 


GND 


Ibata et al. ri997) 






-2.83 ±0.20 


-1.33 ±0.20 


GND 


Dinescu et al. f2005) 


LMC 


50 


±2.03 ±0.08 


±0.44 ±0.05 


HST 


Kalhvavalil et al. f2006b) 


SMC 


57 


±1.16 ±0.18 


-1.17±0.18 


HST 


KaUivavalil et al. f2006a) 


UMi 


68 


±0.5 ±0.8 


±1.2 ±0.5 


GND 


Scholz & Irwin ri994) 






±0.056 ±0.078 


±0.074 ±0.099 


GND 


Schweitzer et al. fl997) 






-0.50 ±0.17 


±0.22 ±0.16 


HST 


Piatek et al. r2005) 






-0.04 ±0.07 


0.13 ±0.07 


MEAN 




Sou 


79 


±0.72 ±0.22 


-0.06 ±0.25 


GND 


Schweitzer et al. fl995) 






±0.09 ±0.13 


±0.02 ±0.13 


HST 


Piatek et al. f2006) 






±0.25 ±0.11 


0.00 ±0.12 


MEAN 




Dra 


82 


±0.6 ±0.4 


±1.1 ±0.5 


GND 


Scholz & Irwin fl994) 


Car 


103 


(±0.22 ±0.09 


±0.15 ±0.09) 


HST 


Piatek et al. f2003)'' 






±0.22 ±0.13 


±0.24 ±0.11 


HST 


Piatek (2007, priv. comm) 


For 


140 


±0.59 ±0.16 


-0.15 ±0.16 


GND 


Dinescu et al. f2004) 






±0.476 ± 0.046 


-0.36 ±0.041 


HST 


Piatek et al. f2007) 






±0.485 ± 0.044 


-0.354 ±0.028 


MEAN 





"Proper motion measurement without advanced CTI correction for HST/STIS. 
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Table 2: Shape parameter 7 and strength parameter C, of the distribution of directions of 
normal vectors fitted to 10 000 bootstrap samples of the simulated CDM satellite galaxies. 
The first column gives an identifier of the sir nulated haloes. Fo r comparison, the values for 
the Milky Way are given in the last row. See iMetz et al.l (120071 ) for further details. 



shape parameter 7 strength parameter ^ 



ghl 
gh2 
gh3 
gh6 

gh7 

ghlO 

MW 



0.5 
2.3 
0.3 
1.2 
3.9 
1.5 
3.85 



3.2 
4.2 
3.5 
3.0 
2.6 
2.8 
4.29 
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Fig. 1. — Orbital poles of Milky Way satellite galaxies as derived from their measured 
proper motion and radial velocities. The directions are shown in an aitoff projection in 
Galactocentric coordinates. The solid lines give the projected arc-segments derived from 
the uncertainties of the measured proper motions. Different symbols mark data derived 
by different methods: circles: HST, diamonds: ground-based measurements, and hexagons: 
the weighted mean values from Table [H The filled star symbol marks the mean spherical 
direction / of the directions of the angular momenta of the satellites excluding Sagittarius 
and Sculptur, and the solid loop gives the spherical standard deviation of this sample. The 
smaller, open star symbol marks the mean spherical direction as before, but now treating the 
LMC/SMC as a bound system whose barycentre is moving with the velocity of the LMC, 
assuming an LMC/SMC mass ratio of 5/1. The dashed loops indicate regions with 15° 
and 30° from the direction of the normal to the plane fitted to the 11 classical Milky Way 
satellites (the DoS pole). Note the proximity of I to the normal of the DoS. 
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Fig. 2. — An aitoff projection of the locations of the orbital poles from the six CDM simu- 
lations. The projection is centred on the potential rotational axis - the normal of the fitted 
plane - if the formally fitted disc of satellites were rotationally supported for the simulated 
data (see text for details). The dashed- loops give the same regions as in Fig. [H the star 
symbol the spherical mean direction of the angular momenta, which is nowhere closer than 
~ 30° to the normal of the fitted planes. 
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Fig. 3. — Cumulative distribution function of the pole distance function determined for 
the Milky Way and the CD M simulated satellite galaxies derived for six MW-halo models 
from iLibeskind et al.l (120051 ) . The continuous lines give the theoretical measured FDD for 
8 (MW), 11 (CDM), and the dotted line shows the case for an infinite number of random 
directions (Eqn. [2]). 



